Induced interactions in dilute atomic gases and liquid helium mixtures 
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In dilute mixtures of two atomic gases, interactions between two minority atoms acquire a con- 
tribution due to interaction with the majority component. Using thermodynamic arguments, we 
derive expressions for this induced interaction for both fermions and bosons for arbitrary strength 
of the interaction between the two components. Implications of the work for the theory of dilute 
solutions of ''He in liquid *He are discussed. 
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I. INTRODUCTION 

Induced interactions are responsible for a variety of 
phenomena in condensed matter physics, ranging from 
the superconductivity of metals, to the stability of the 
A phase of superfluid liquid ^He, to the effective inter- 
actions between ^He atoms in liquid '^He [T]. Recently, 
following the work of Mora and Chevy f51 , it was shown 
that such processes are important in dilute mixtures of 
atomic gases [3]. In this paper, we expand on Ref. [3], 
which considered two fermonic species, and extend the 
results to boson-fermion and boson-boson mixtures. In 
addition, we draw a number of conclusions relevant to 
dilute solutions of ^He in liquid ^He. 



II. FERMION-FERMION MIXTURES 

We begin by giving a compact derivation of the basic 
result that the Landau quasiparticle interaction for two 
minority atoms (denoted by 2) in a Fermi gas consisting 
of a majority species 1, has a momentum independent 
contribution of the form 
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is the number of majority atoms that must be added per 
minority atom in order to keep the chemical potential 
of the majority species fixed and dni/djii = A/i(0) — 
TOiPi?i/27r^ is the density of single-particle states at the 
Fermi surface for species 1. (We use units in which h is 
equal to unity.) The density of species i is denoted by n^, 
the mass of an atom by mi , and the Fermi momentum by 
PFi- Physically, the quantity v is the number of majority 
atoms in the dressing cloud of a minority atom. 



To derive the result ([T]) , it is simplest to work in terms 
of an effective low energy theory, in which high-lying 
states are eliminated, and only low-lying ones are re- 
tained. The effective low-energy Hamiltonian is 
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where V is the volume of the system, iVi is the total num- 
ber of 1-atoms, £2 is the energy to add a single 2-atom to 
the 1-atoms, is the effective mass of a single 2-atom 
and g is the strength of the effective interaction between 
different atoms. The prime on the sum indicates that 
the q = term is omitted, since this is included in the 
£2 term. We shall assume that the momentum scale for 
variations of the coupling strength are large compared 
with the Fermi momentum, and therefore the momen- 
tum dependence of the coupling may be neglected. Also, 
we have not written explicitly the direct interaction be- 
tween two 2-atoms or between two 1-atoms because this is 
short-ranged and, consequently, the direct and exchange 
contributions to the interaction energy cancel. We shall 
consider only the isotropic part of the interaction and ne- 
glect the momentum dependence of the interaction which 
leads to contributions to Landau parameters other then 
that for / = 0. 

The coupling constant g gives the change in the in- 
teraction energy between small long-wavelength density 
disturbances in the medium. In the limit of small concen- 
trations of 2-atoms the contribution to the energy from 
the Fermi motion of the 2-atoms is negligible and there- 
fore the interaction energy may be replaced by the total 
energy. Thus one sees that [1] 
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where £{ni,n2) is the energy density as a function of the 
densities of the two components. 

The induced interaction is always at least of second 



order in g and we evaluate it by calculating the con- 
tribution to the energy, which is given by 
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where /p is the particle distribution function, the index 
on the momentum variable indicating whether it refers 
to 1-atoms or 2-atoms. From Eq. ^ one may calculate 
the corresponding contribution to the Landau effective 
interaction between two 2-quasiparticles, 
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where P2 and P2 arc taken to be vanishingly small. Thus 
one finds 
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This is positive, since although the contribution to the 
total energy is negative, one of the distribution func- 
tions corresponds to a hole line, which carries a factor 

1-/] 
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and therefore the second functional derivative 
of the energy with respect to the distribution function for 
2-atoms is positive. Expressed in the language of Ref. [5], 
Pauli blocking reduces the magnitude of the negative con- 
tribution to the energy, thereby giving a positive contri- 
bution to /. An equivalent description is that, while the 
induced interaction is intrinsically negative, its contribu- 
tion to the Landau quasiparticle interaction comes from 
an exchange term, which gives an additional minus sign 
[3]. Contributions of higher order in g will contain ad- 
ditional powers of the Fermi momentum of the minority 
component and are therefore negligible in the limit of a 
low concentration of minority atoms. 



III. A BOSE GAS WITH A DILUTE FERMI 
COMPONENT 

An analysis similar to that in Sec. |IT] may be carried 
through for a majority Bose component. The elementary 
excitations in the Bose system are phonons and the low- 
energy effective Hamiltonian is 

H =E„{Ni) + sqa^a^ + e2A^2 + ^ 
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where the operator ajj creates phonons in the Bose sys- 
tem and Eo{Ni) is the energy of A^i bosons includ- 
ing the effects of interactions. Again we neglect the 
momentum-dependent part of the effective interactions 
due to, e.g., couphng of fermions to the superfluid ve- 
locity of the bosons. The energy of a low-wavelength 
phonon of wavenumber q in the Bose system is sq, where 
= {ni/mi)d^i/dni. The matrix element il/g for the 
density operator to create or destroy a phonon of mo- 
mentum q is given by [5] 



mq_ 
2mis 



1/2 



(9) 



There are no terms with higher powers of a and since 
phonons in a Bose-Einstein condensed gas exhaust the 
frequency-weighted sum rule for the density-density cor- 
relation function. The contribution to the energy calcu- 
lated in second-order perturbation theory is 
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where riq is the phonon distribution function. Thus for 
two fermions on the Fermi surface, the effective interac- 
tion given by Eq. ^ is 
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This result is valid regardless of the strength of the boson- 
boson interaction. The only difference compared with the 
case of a majority gas of fermions is that Mi ~ dni/d^i 
in Eq. ([7| must be replaced by the expression for a Bose 
gas. For a weakly-interacting Bose gas, /ii = niUn, 
where Un = Anan/mi is the effective low-energy in- 
teraction, flu being the scattering length, and therefore 
dfii/dni = Uii. 



IV. MINORITY BOSE COMPONENT 

Analogous arguments may be carried through for a 
mixture with a minority Bose component. We consider 
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the case when the bosons are in a Bose-Einstein con- 
densate and do not treat the case where bosons form 
fermionic diatomic molecules with the majority fermions 
[51 [7] . The bosons may thus be described by their density, 
n2, and the boson superfluid velocity, v^. For a spatially 
uniform system consisting of bosons, its ground state en- 
ergy is a function only of the boson density, n2, so the 
effective interaction is defined by the relation 
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For moving condensates, there will also be effective in- 
teractions involving the superfluid velocity, but we shall 
not take these into account explicitly here. One differ- 
ence compared with the case of minority fermions is that 
the direct boson-boson interaction does not vanish. As 



defined in Eq. (12 1, there is no exchange term since the 



2-atoms are in a Bose condensate. Consequently the ex- 
change process is identical to the process without ex- 
change and to include it explicitly would be double count- 
ing. A second difference is that the induced interaction 
contribution to / is negative, because the wave function 
for bosons is symmetric under interchange of particles 
and, consequently, the extra minus sign acquired in the 
case of fermions is absent. The result is 
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The collective modes of the Bose system may be cal- 
culated by the Bogoliubov approach in which the effec- 
tive interaction between bosons is given by /, provided 
the frequency dependence of the interaction may be ne- 
glected. This condition is satisfied provided the velocity 
of the Bogoliubov mode is much less than the Fermi ve- 
locity or the sound speed of the majority component. 
Thus the sound speed, S2, in the Bose gas is given by 
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If / is negative, the sound speed is imaginary and density 
modes in the system are unstable. This corresponds to 
the thermodynamic stability condition, see, e.g., Ref. ^8; 
Sec. 12.1.1] 
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since /^"' = d^2/dn2. 



V. DYNAMICAL EFFECTS 

Here we show how density modes are affected by the 
induced interaction. For definiteness, let us consider the 
response of the density of a minority fermion component 
to a potential acting on the minority atoms. We assume 
the wave vector q of the perturbation to be small. The 
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FIG. 1: Diagrams representing the induced contribution to 
the effective interaction / (left) and the corresponding dia- 
gram for the crossed channel (right). The wavy line repre- 
sents exchange of an excitation in the majority component, 
either a phonon in the case of a Bose gas or a particle-hole 
pair in a Fermi gas. 



contributions to the effective interaction between minor- 
ity atoms have the forms shown in Fig. ([T]) : the left-hand 
diagram corresponds to the contribution to the effective 
interaction that we have calculated earlier, while the sec- 
ond term is due to the response of the majority atoms 
at wave vector q. Summing up all bubble diagrams, the 
density-density response function for the minority atoms 
is given by 
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/(q,^)=52^i(q,w). (18) 
If the majority atoms are fermions, xi{^.i^) is given by 
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while if they are bosons it is given by 
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In the limit w/g — )■ 0, /(q, w) tends to /. The can- 
cellation of the two terms is a consequence of the anti- 
symmetry under interchange of two fermions in the same 
internal state and the fact that for uj/q Q the contri- 
butions of the two processes shown in Fig. [T] are equal in 
magnitude. Consequently, effects of the 1-2 interaction 
disappear in the density response. However, at nonzero 
frequency, there will in general be effects due to the in- 
teraction because the ratio of the energy transfer to the 
momentum transfer is different in the two particle-hole 
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channels for the 2-atoms. In words, while the Pauli prin- 
ciple forbids two atoms in the same internal state being 
at the same point in space at a given time, it does not 
forbid two such particles being at the same point at dif- 
ferent times. 

The frequency dependence of the interaction has strik- 
ing implications for the scattering rate of 2-atoms by 2- 
atoms, since for a; = the effective interaction vanishes 
and consequently the scattering rate will grow with tem- 
perature as T^, in contrast to the behavior predicted 
by standard Landau theory f9i . Related effects occur for 
quark-gluon plasmas, where again the frequency depen- 
dence of the scattering amplitude plays a decisive role 
in determining scattering rates at low temperature [10]. 
However, the rate for scattering of 2-atoms by majority 
fermonic atoms will still have the usual dependence. 

For a minority Bose gas, the frequency dependence of 
the effective interaction affects the dispersion relation of 
the Bogoliubov mode, which is given by 
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For a majority Bose gas, x(q, w) increases for small uj 
and consequently the sound speed is less than the result 
( 14 1 predicted on the basis of the static interaction ( 13 ). 



The situation for a majority Fermi gas is different, and 
for small uj/viq, 
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The leading effect of dynamics at low frequencies is an 
imaginary contribution to xi due to Landau damping. 
Thus sound waves in the Bose gas can decay into particle- 
hole pairs of the majority species. The time t for decay 
of the intensity of the wave is given by 
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which is small compared with the real part of the fre- 
quency for small concentrations of the minority compo- 
nent. In Appendix A we present on the basis of a func- 
tional integral approach an alternative derivation of a 
number of results in this section for the case of dilute 
gases. 



the potential are typically obtained by assuming a par- 
ticular form for the spatial Fourier transform of the inter- 
action and then fitting parameters to obtain agreement 
with measured transport coefficients. In this approach, 
the dynamics of the ^He atoms is taken into account im- 
plicitly, since it is assumed that the ^He atoms respond 
on a time scale short compared with characteristic times 
for the '^He. In the Landau theory of the dilute solutions, 
the response of the ""^He is taken into account through its 
effect on contributions to the Landau parameters. This 
effective interaction corresponds to the definition 
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where /p is the distribution function for '^He quasipar- 
ticles and /i4 is the ^He chemical potential. With this 
definition, the local density of ^He atoms adjusts to the 
local density of "^He atoms, and it is this interaction that 
should be used in formulating a theory of the static prop- 
erties of mixtures, such as the magnetic susceptibility or 
the ■^He contribution to the bulk modulus, if one wishes 
to avoid treating explicitly the response of the '^He. 

For calculating dynamical effects, the "^He density must 
to be treated as a dynamical variable, and therefore the 
appropriate definition of an effective interaction is 
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This corresponds to the natural generalization of the orig- 
inal Landau definition of a quasiparticle interaction to 
allow for a second component. For the case of a major- 
ity Fermi gas rather than "^He, the density and the 
superfiuid velocity V4 would be replaced by the distri- 
bution function for the majority component. Dynamical 
processes may then be described in terms of, e.g., the 
kinetic equation for "^He quasiparticles and the equations 
of superfiuid hydrodynamics for the ^He. 

We now consider the relationship of the two different 
definitions of effective interaction to microscopic theory 
[TT] . The standard definition of the quasiparticle inter- 
action in a normal Fermi system is |12j 

/p,p. = lim lim z^T{p,p\q,uj), (26) 



VI. IMPLICATIONS FOR DILUTE SOLUTIONS 
OF HELIUM ISOTOPES 

Our considerations above have implications for the the- 
ory of dilute solutions of '^He in liquid ^He. The standard 
approach adopted by Bardeen, Baym and Pines [T](BBP) 
is to assume that the interaction between two ^He impu- 
rities may be modeled by a potential that is local in time 
but of nonzero range in space [Tj . This potential includes 
effects of the induced interaction between ■^He atoms due 
to exchange of excitations in the ''He. The parameters of 



where z is the wave function renormalization parameter, 
and r(p,p',q, cj) is the two-particle vertex function for 
bare particles. The momenta p and p' are those of the 
incoming particles, uf is the Fermi velocity, q is the mo- 
mentum transfer and lo the energy transfer. In the case 
of mixtures, the vertex function depends not only on the 
distribution of fermions but also has contributions from 
interactions of fermions with bosons. If one applies the 
definition (26 1 to helium mixtures, contributions to the 
vertex function for scattering of two '^He atoms due to 
exchange of phonons in the Bose system carrying mo- 
mentum q will vanish, because the phonon propagator 
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behaves as n/^q^ jm^ij^ (cf. Eq. (20)). In physical terms, 
the ^He does not respond because of the high frequency. 
For a low concentration of ^He, vp is very much less 
than s. The effective interaction used by BBP includes 
the effects of phonon exchange in the particle-hole chan- 
nel with momentum q and in the microscopic theory it 
corresponds to the quantity 



/, 



BBP 
Pt.p'i 
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In physical terms, w/g must be much less than s in or- 
der that the response of the '^He atoms be given by its 
zero-frequency value. In the BBP approach, exchange is 
taken into account explicitly, and therefore, in defining 
an effective interaction, the two "^Hc atoms are taken to 
be in different spin states. 

In the standard definition of effective interactions for 
a Fermi liquid, the two particle-hole channels are treated 
differently in the case of two fermions in the same spin 
state. For the channel with momentum transfer q, the 
limit taken is the high frequency one, while for the chan- 
nel with momentum transfer p + q — p', it is the low 
frequency one, since the energy transfer is zero. The 
effective interaction vanishes for p = p' in the Born ap- 
proximation but not when many-body processes involv- 
ing particle-hole pairs are taken into account. This prob- 
lem has been investigated in detail in Ref. [13]. 

VII. CONCLUDING REMARKS 

In this paper we have shown how the contribution 
to the effective interactions between minority atoms in- 
duced by interactions with a majority component may 
be expressed in terms of thermodynamic quantities. An 
important remark is that the results are not limited to 
cases where the majority gas is weakly interacting. Ex- 
perimentally, an interesting case to explore would be a 
Bose gas with weak, repulsive interactions. In this case 
the sound speed is small and therefore the effects of the 
induced interaction can be correspondingly large. 

In addition to the interaction mediated by coupling 
to density fluctuations in the majority component, there 
is also an interaction between minority atoms due to 
coupling via a current-current coupling [1] or, in the 
case of a fermion majority component, distortions of 
the majority Fermi surface that are not spherical. In 
general these too should be included, but they are 
generally much smaller than the interaction induced 
by the density-density coupling. This interaction does 
not contribute to the effective interaction between 
two fermions at the Fermi surface in a majority Bose 
component because the current operator for transitions 
of a fermion between two states on the Fermi surface is 
transverse, and therefore it cannot create a phonon in 
the Bose gas, which is longitudinal. 
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Appendix A: Functional integral approach to 
Fermi-Bose mixtures 



Here we give a derivation of some of the results in 
Sec. V based on a functional integral approach. We con- 
sider a binary mixture consisting of a fermion species and 
a boson species at zero temperature and we shall assume 
that the densities of the two components are so low that 
binary interactions dominate. The Hamiltonian is 
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where the field operators are </) for the bosons and "i/; for 
the fermions and the integral is over coordinate space. 
The interactions between bosons are repulsive. The cou- 
plings are given by qbb = ^TraBB/mB > and jij 



l/imrgBp) + A/tt^ = l/2TTaBF- 
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Here the reduced mass is l/m^ — l/rriB + 1/mp, A is 
the momentum cutoff, and ubb a-nd gbf a-re the s-wave 
scattering lengths. 

To simplify the discussion, we assume that the Bose 
gas is dilute, in the sense that a^^ns ^1. In the 
case that the fermions are the minority and the Fermi 
momentum kp ^ 0, the induced interactions between 
the minority fermions mediated by the majority bosons 
gives rise to a nonzero Landau parameter as discussed 
in Sec. III. In the case that the bosons are the minor- 
ity, we shall consider situations where bound states of 
bosons and fermions (fermionic diatomic molecules) need 
not be taken into account. This could be due to either the 
boson-fermion interaction being so weak that there are 
no bound states or to the state under consideration being 
a metastable one in which there are no such molecules. 
When bound states are present, the perturbation theory 
treatment given here needs to be extended to allow for 
the presence of molecules. In the path integral represen- 
tation, after the fermions been integrated out and terms 
up to the second order of qbf retained, the partition 
function is given by 
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Here the effective action for the bosons is given by 

72 
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where the symbol 1 stands for {ri,ti}. The integral 
over imaginary time t is from to — i/3, and — 
{2mpiip)^^^ /Qn^ . The induced interaction has the form 
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where G^'* is the free fermion Green function. 

We apply the Bogoliubov approximation for the bosons 
and obtain for the frequency Q of bosonic modes the 
dispersion relation 
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where P is the momentum of the mode and ^ind is the 



Fourier transform of Eq. ( A5 1 . In the long wavelength 
limit P/kp 0, 
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Here vp = kp/mp is the Fermi velocity. The sound speed 
c = n/P for P — >■ and — >■ + iJ and in the regime 
c/vp <^ 1 is 
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with cq = \/ gBBnB /mB- The reduction of the sound 
speed, Rec < cq, is due to the attractive induced inter- 
actions between bosons mediated by fermions [TH [T5] . 
The small imaginary part of c indicates that the modes 
decay into particle-hole pairs of the majority fermions 
(Landau damping). 

Equation ( |A8[ ) implies dynamic instability for 

gBB ~ gBP'^^F /^f^P — 0- The same instability con- 
dition can be deduced from the energy density of the 
mixture 
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where the first term is the energy of a filled Fermi sea, by 
requiring that the variation of £ to second order in the 
density variations be negative [5]. Note that Eq. (AlO) 
takes into account the interaction energy only at the 
mean field level. 

When gsp is no longer small, a similar dynamic insta- 
bility condition can be derived following the above argu- 
ment. We can still formally integrate out the fermions 
in the partition function and obtain an effective action 
for the bosons. Since hb — > 0, in the effective action 
for the bosons, the interaction effects due to gBF can 
be taken into account in two steps. The first step is to 
modify the properties of a single boson, such as the effec- 
tive mass m*g, the single boson energy fip (as for polar- 
ized fermions) and the quasiparticle residue z{< 1). The 
second step is to change the interactions between these 
bosonic quasiparticles. In the low energy and long wave- 
length limit, induced interactions between the bosons are 
given by the diagrams shown in Fig. ([T]) with the vertex 
replaced by d^p/ drip'' . Within the Bogoliubov approxi- 
mation, we conclude that if 
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the system becomes dynamically unstable. 
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